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Abstract 

A nonlinear kinetic equation for nonrelativistic quantum plasma with electromagnetic interac- 
tion of particles is obtained in the Hartree's mean-field approximation. It is cast in a convenient 
form of Vlasov-Boltzmann-type equation with "quantum interference integral" , that allows for rela- 
tively straightforward modification of existing classical Vlasov codes to incorporate quantum effects 
(quantum statistics and quantum interference of overlapping particles wave functions), without 
changing the bulk of the codes. Such modification (upgrade) of existing Vlasov codes may provide 
a direct and effective path to numerical simulations of nonlinear electrostatic and electromagnetic 
phenomena in quantum plasmas, especially of processes where kinetic effects are important (e.g., 
modulational interactions and stimulated scattering phenomena involving plasma modes at short 
wavelengths or high-order kinetic modes, dynamical screening and interaction of charges in quan- 
tum plasma, etc.) Moreover, numerical approaches involving such modified Vlasov codes would 
provide a useful basis for theoretical analyses of quantum plasmas, as quantum and classical effects 
can be easily separated there. 
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I. INTRODUCTION 



The interest in quantum plasmas - plasmas with quantum effects playing a significant 
role in their collective behavior - has considerably increased in the recent decade, during 
which a significant number of publications appeared on this subject (see, e.g., lMl2| and ref- 
erences therein). This surge of interest can primarily be associated with the recent progress 
in manufacturing and manipulation of metallic and semiconductor nanostructures, whose 
properties are to a large extent governed by collective (plasma) effects of their electron 
(and hole) population. Quantum plasma effects may also become important in fast ignition 
scenario of modern inertial confinement fusion (ICF) experiments, when deuterium-tritium 
plasma is compressed to very high (super-solid) densities b y in tense laser pulses 13Ml5| . Re- 
cent spectral measurements of x-ray Thomson scattering 



electron distribution function in warm dense matter regime 18| , which opens the possibility 



17( enable precise probing of 



of direct experimental studies of weakly and strongly coupled quantum plasmas. All these 
developments make quantum plasma regime both accessible and relevant, making theoreti- 
cal efforts in describing quantum plasmas, and in particular linear and nonlinear collective 
processes in quantum plasmas, both timely and important. 

Microscopic theoretical models of quantum plasma usually treat the plasma as a weakly 
coupled ga s of particles, described in the framework of the Hartree's mean-field approxi- 



mation 



19|. We should note that not all quantum plasmas are weakly coupled, e.g., for 



conductivity electrons in metals the potential energy of their interaction with each other 
and with the lattice ions is comparable to or exceeds their mean kinetic energy, thus mak- 
ing them a moderately coupled quantum liquid for which the mean-field approximation is, 
strictly speaking, invalid. However, even for such quantum plasmas the simple mean-field 
collisionless approximation yields a good qualitative picture of many collective phenomena 
such as waves and their interactions, while accounting for particle correlations (e.g., due to 
collisions and exchange interaction of particles) only yields quantitative corrections to the 
results of the mean-field theory; see, e.g., Ref. jsO | - For this reason, here for simplicity we 
restrict to the mean-field collisionless approximation for quantum plasmas. This assump- 
tion is quantitatively justified for plasmas with a small coupling parameter Q V <C 1, and 
captures most qualitative features of collective modes even for plasmas with T ~ 1, where 
r = U int /e kin ~ e 2 n 1 / 3 /ekin, ^int an d ekin are the mean potential and kinetic energies of 
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plasma particles, e is the particle charge, n is the particle number density. 



The linear theory of weakly coup 
tized, has been developed by 1960s 



ed 
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quantum plasma, both unmagnetized and magne- 



23} . However, quantum plasmas are often formed 



or exist under the conditions of powerful energy input (e.g., in ICF experiments), leading 
to a rapid development of nonlinear phenomena, which require a proper theoretical treat- 



ment. A number of papers (see, e.g., Refs [2|, [ll|, l24j427j| and references therein) appeared 
recently addressing various nonlinear phenomena in quantum plasmas by employing the so- 
called quantum fluid theory (QFT) approach [28 31 ] . The QFT is considerably simpler 
than the more general kinetic theory (e.g., the theory based on the kinetic equation for the 
quantum plasma distribution function, the so-called Wigner function 321]). as QFT operates 
with quantities (fluid density, flow velocity, pressure) that depend on four variables r, t, 
while the kinetic theory operates with the quantum distribution function / that depends 
on seven variables r, p, t, where r and p are position and momentum vectors, respectively, 
and t is time. This relative simplicity makes the QFT description more favorable (compared 
to the kinetic description) when treating nonlinear phenomena, however, it comes at a cost 
of rather restricted validity range of such description and, correspondingly, its results (the 
characteristic length scale of waves or nonlinear structures should be large compared to 
the Thomas-Fermi length in degenerate quantum plasma, or to the Debye length in non- 
degenerate quantum plasma, for QFT to be valid). This restriction appears as a result of the 
assumptions necessarily made in the course of derivation of the "collisionless QFT" from the 
kinetic theory of collisionless quantum plasma jsl, [^sj. Also, as all fluid theories, QFT does 
not account for essentially kinetic effects, such as, e.g., Landau damping [33|. Because of its 
restricted validity, QFT is unable to correctly describe nonlinear wave interaction processes 
such as, for example, stimulated Raman scattering (SRS) and stimulated Brillouin scattering 
(SBS) of laser light off electrostatic plasma modes at short wavelengths (note that, unlike in 
classical plasmas, in degenerate quantum plasma electrostatic modes can have zero damping 
at short wavelengths |8(, and thus can effectively participate in nonlinear interactions of 
electrostatic and electromagnetic waves), or other modulational-type interactions involving 
higher-order kinetic modes in quantum plasma. Hence, a kinetic model is needed for proper 
description of such processes. Proper description of dynamical screening and interaction of 



charged bodies in quantum plasma also requires a kinetic treatment 



fly 



35|. 



In this paper, we develop, in the framework of Hartree's mean-field approximation 19[, a 
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kinetic model of a nonrelativistic collisionless quantum plasma consisting of spinless charged 
particles with electromagnetic interaction. This model covers a wider range of plasma pa- 
rameters than the classical kinetic model based on Vlasov equation coupled with Maxwell's 
equations, and should be useful for studying a wide range of linear and nonlinear electro- 
static and electromagnetic processes in unmagnetized or magnetized uniform or nonuniform 
quantum plasmas, for which the effect of collisions, particle spin and relativistic effects are 
unimportant or only lead to minor quantitative corrections. (Hence this model is not ex- 
pected to be applicable to, e.g., SRS of relativistically strong laser pulses in plasma, and 
is not capable of describing, e.g., spin waves). The main result is the kinetic equation for 
the one-particle quantum distribution function (one-particle Wigner's function fx) of such 
system, which is presented in the form of a "classical" Vlasov-Boltzmann equation for j\ 
(which by itself contains quantum effects), with explicit quantum terms appearing on the 
right-hand side in a "quantum interference integral" (manifesting the effect of quantum in- 
terference of overlapping particle wave functions, and vanishing in the classical limit). This 
form of the quantum kinetic equation allows for its rather straightforward implementation 
in existing numerical Vlasov codes, whose only modifications would be the inclusion of the 
"quantum interference integral" as a source term in the Vlasov equation, coupled with the 
correct initial and boundary conditions for the quantum distribution function (the latter in 



32|]). Exist- 



general can be different from those for the classical distribution function, see 
ing classical Vlasov codes, modified in this way, can be used to simulate various linear and 
nonlinear collective effects in nonrelativistic quantum plasmas. 



II. KINETIC MODEL OF QUANTUM PLASMA WITH ELECTROMAGNETIC 
INTERACTION 

In this section, we start with recalling the well-known basics of the quantum statistical 
theory, and then proceed to obtain the kinetic equation for a quantum plasma with spinless 
particles interacting electromagnetically, which is the main result of this paper. 
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A. Quantum distribution function 



Statistical description of a system of many interacting quantum particles (quantum 
plasma) is done in terms of its density matrix, which allows us to obtain the mean val- 
ues and probability distributions of all the physical parameters of the system. In particular, 
it is convenient to describe the quantum plasma in terms of the density matrix in the mixed 
coordinate-momentum representation - the quantum distribution function /jv(rjV; Pjv, t), 
suggested by Wigner [36( and sometimes called the iV-particle Wigner function, defined as 

1 



/at (r n, PN,t) 



(2tt) 



3A ! 



dX e 



-i\-p N 



Pn(*n - h\/2, r N + hX/2, t), 



(1) 



where pN(r^,r' N ,t) is the density matrix of the system in the coordinate representation, N 
is the number of particles in the system, r^v and pn are 3 N- dimensional vectors denoting 
the sets of coordinates and canonical momenta of all system particles, and h is the reduced 
Planck constant. In the classical limit h — > 0, Jjv becomes the classical iV-particle phase 
space distribution function, hence the description in terms of the Wigner function covers 
both classical and quantum plasmas. The properties of the Wigner function are discussed in 



detail by, e.g., Tatarskii 



32|. We note that the Wigner function (JTJ) should not be thought 



of as the density of the system's states in r^r, Pn phase space, as it is not positively definite 
(i.e., can attain negative values) due to non-commutativity of coordinate and momentum 
operators in quantum mechanics 32] . 

Equation governing the evolution of /jv is obtained from the evolution equation for the 



density matrix pn(tn, r'^, t) and reads 
df N (r N ,p N ,t) 1 



i9i,y. 



3 



at 



(2ny N h 

x 



dX N dk N dff N dq N e 4^-^-P-)+ k --(^- r -)]^( q7V , ff N , t) 



N 



q N + ^hX N ,ff N - -Hk N 



(2) 



where ^(r^r, p^v, t) is the system's Hamiltonian function containing exact (not averaged) 
fields through which particles interact. 

Description of plasma in terms of the iV-particle quantum distribution function 
f]v(r N ,p N ,t) is very cumbersome, as /jy depends on a huge number of variables, 6iV + 1. 
However, for most physical applications it is sufficient to know the one-particle quantum 
distribution function /i(ri,pi,t) that depends on coordinates ri and momenta pi of one 
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particle, irrespective of positions and momenta of all other particles in the system: 
/i(ri,pi,t)= / f N (r N ,p N ,t)dr 2 . . .dr N dp 2 . . .dp N . 



The equation for f\ can be obtained from ([2]) using the approach similar to the BBGKY 
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(Bogoliubov-Born-Green-Kirkwood-Yvon) hierarchy approach for classical plasma 
the resulting equation for pi, t) contains the two-particle Wigner distribution function 
/2(ri, r 2 , Pi, P2j t), the equation for which in turn contains the three-particle distribution 
function, and so on. The resulting set of coupled equations - a quantum analogue of the 
BBGKY hierarchy - is equivalent to Eq. fl2]), and its solution is equally difficult. How- 
ever, in systems of weakly interacting particles, with T = f/intAint <S 1, the problem of 
finding f\ can be significantly simplified by neglecting the two-particle correlation function 



02 ( r i) r 2, Pi 5 P2, t) defined by the following relation [l9j: 



/ 2 (ri,r 2 ,Pi,P2,t) = A(ri,pi,t)/i(r 2 ,p 2 ,t) +^ 2 (ri,r 2 ,pi,P2,t). (3) 

The function g 2 characterizes statistical correlation between particles 1 and 2 due to their 
interaction. It defines the value of the collision integral I(ti, pi,t) in the evolution equation 
for fx, and hence neglecting g 2 in the case of weak particle interactions corresponds to the 
collisionless plasma approximation, J(ri,pi,t) = 0. The resulting approximate equation 
for fi(ri,pi,t) does not contain f 2 , having the form of Eq. ([2]) with N — 1, in which the 
Hamiltonian function K(r 1; pi, t) now contains the mean self-consistent fields averaged over 



mean-field approximation 



3. 



fi, instead of the exact unaveraged fields 19]. Such approximation is called the Hartree's 



B. Kinetic equation for one-particle quantum distribution function 

We consider a gas of spinless charged quantum particles with charge e and mass m (e.g., 
electrons with neglected spin) , whose net charge is compensated by a background of heavy 
immobile particles (e.g., ions) with density no- Neglecting two-particle correlations due 
to exchange interactions and collisions between particles (see the discussion of collisionless 
mean- free approximation in Introduction), the quantum kinetic equation for the one-particle 
Wigner function fi of plasma particles is obtained from ([2]) for N = 1 and has the form (in 
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what follows, we drop the index 1 of the one-particle Wigner function, for brevity) 
dJ ^T 1 = (2^)4 /•••/ ^^qe < P^ P ) +k ^]/(q^,t) 

x nL-hx,ff+hk,t^ -nU+^h\,7f-~hk,tj , (4) 

U(r,P,t) = -L(p-^A(r,t)) 2 + e0(r,t), (5) 



where 4>(r,t) and A(r, t) are respectively scalar and vector potentials of the self-consistent 
electromagnetic field. Equations for <f) and A are obtained from the Maxwell's equations, 
and their form depends on the chosen gauge condition. For example, for the Coulomb gauge 
V ■ A = 0, we have the following coupled equations for <fi and A: 

- V 2 = 4vr[en e (r,t) - en ] , (6) 
1 d 2 A ldV(j) 4tt. , N 

where the electron charge and current densities are defined as the zeroth- and first-order 
moments of the Wigner function: 



en e (r,t) = e J f(r,P,t)dP, (8) 

j e (r,t) = - f p/(r,P,*)dP. (9) 
m J 

Here P = p+(e/c)A(r, t) is the canonical momentum, and p = mv is the kinetic momentum 
of a quantum particle (we use CGS units). Equations (SJ) ([9]) form a closed set of equations 
describing electrodynamics of quantum plasma with electromagnetic interaction of particles. 

C. Reduction to the Vlasov-Boltzmann-type kinetic equation 

Eq. 01]) with (jSJ) can be cast in the form of the Vlasov-Boltzmann-type kinetic equation 
in which all the explicit quantum terms are grouped separately, by the following steps: 



Change of variables in (|I])-(J5]) from r, P, t to r, p, t, noting that 

df(r,P,t) = df(r,p,t) _ e dA(r,t) <9/(r,p,t) 
dt dt c dt dp 



(10) 



The next step is to change the integration variable ff to £ = if — (e/c)A, and perform 
the integration over q and k in the right hand side of (jl]). One thus obtains the 



following intermediate equation for /(r, p, t): 



df f p e ^\ df e dA df e 
dt Vm mc J dr c dt dp mc 



( P + lA 



d_ 
Or 



A 



df_ 
dp 



1 1 

(2vr) 3 m 



<A(C-p) 



-iem/(r,£,t)4 x [0(r,t)] 



c/Ac/£ e 



+-/(r, e, *) (A + - [2A(r, *) - A(r_, t) - A(r+, • D r>X [A(r, t)] 



+ ^ (V/(r, t) + /(r, £, t) Vj • [A(r_, t) + A(r + , *)] 
[([A(r_, f) + A(r + , f)] ■ V) A(r, t)] • 



(11) 



where the "spatial difference" operator D r ^ is defined as D ?[g(r,t)] = 
[g(r + ,t) — g(T-,t)]/h, with r± = r ± HX/2 and g(r,t) being a scalar or a vector func- 
tion. 



To Eq. ( II II) . we add the following identities: 
1 1 



(2vr) 3 m 



d\d£ e 



iA«-p) 



J(r,£,t)(A-V)(£-A) 



— [p x V x A + (p • V) A] • ft^'P'*) + _L( V • A)/(r, p, t), (12) 
mc op mc 



(27T) 



| f d\d£<f^ {zem(r • V0)/(r,£t)} = -eV0 ■ 



5/ 
<9p' 



(13) 



WT^SS dXd ^~ P) {'a ( v/(r '^ } + /(r '^ )v ) ■ A 

= ■ V/(r, p,t)-—(V- A)/(r, p, t), 

mc mc 



(14) 



1 1 

(27r) 3 m 



<iAd£ e 



<A(€-p) 



?c 



mc 



/(r,£t)A-[(A-V)A] 



[(A • V)A] 



df(r,p,t) 
dp 



(15) 



Collecting terms, and noting that 



d/(r,£,t) 



s iA-(€-p) 



5 



^ .9^ L 

we finally arrive at the desired result shown below. 



iA-(S-p) 
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D. Result: Vlasov-Boltzmann-type kinetic equation with "quantum interference 
integral" 

The above procedure yields the following resulting form of the kinetic equation (JJ|) for 
the quantum distribution function /(r, p, t): 



df , P df 
at m or 



^ pxB 
E+ 



mc 



df 

^ = / ff (r,p,t), (16) 



where E = — V0 — (l/c)dA/dt and B = V x A are the electric and magnetic field strengths, 
respectively. The left hand side of (fT6j) formally coincides with that of the classical Vlasov- 
Boltzmann equation (but for the quantum distribution function, which by itself contains 
quantum effects), and the right hand side respresents the "quantum interference integral" 
I q (r, p, t) defined as 

I q {r,P,t) = 777^ // d\d£exp[i\ ■ (£- p)] 



(2 



7T 



Or 



x { ie f(r,U)[\-^ 1 -D rX [4 > (r,t)] 



2mc 



■[2A(r,t)-A(r_,t)-A(r + ,t)]}. (17) 

The "quantum interference integral" I q manifests the effect of quantum interference of over- 
lapping particle wave functions. It vanishes in the classical limit, which can be seen by 
formally taking the limit h — > in (fT7|) . In this case, the kinetic equation ( 1T5"|) reduces to 
the classical Vlasov equation for the classical distribution function. Note that the kinetic 
equation (fT6|) and the "quantum interference integral" (fT7|) are derived without assuming 
any gauge fixing for the electromagnetic field potentials, and thus can be used in their 
present form for any gauge. 
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III. USE OF THE MODEL 



A. Numerical simulations of nonlinear electromagnetic processes in quantum plas- 



mas 



Quantum kinetic equation (1161) is presented in a convenient form that allows to in- 
corporate quantum effects into existing classical electromagnetic Vlasov (Vlasov-Maxwell) 



codes 



42 



44J, that can then be used for numerical simulation of nonlinear electrodynamic 



processes (e.g., modulational interaction of electrostatic and electromagnetic waves) in quan- 
tum plasma structures. Indeed, the left hand side of (I16p has a form of classical Vlasov 
equation for the quantum distribution function / (which is however not equivalent to the 
classical distribution function), while all the quantum interference effects due to overlap- 
ping of plasma particle wave functions are contained in the right hand side of (1161) . The 
six- dimensional "quantum interference integral" I q {r, p, t) ( JT71) conserves the spatial number 
density of particles: 

J J g (r,p,t)dp = 0, (18) 

which can readily be seen by integrating (|T7|) over p. In the classical limit h — > 0, I q (r, p, t) 
tends to zero as Ti 2 : 

h 2 l 



-Jg(r,P,t) 



24 (2ny 



d\d^ e 



<a-«-p) 



-ief(r, £, t) (A • V)^> + — /(r, e, t) U • (A ■ V)^A 

mc 



3c 



mc 



V/(r,£,t) + /(r,£,t) 



V + — A x (V x A) 



(A ■ V) 2 A \ . (19) 



To include quantum effects into an existing classical Vlasov-Maxwell numerical simulation 
code, only the following two modifications of the code are needed: 

1. Initial equilibrium distribution 

Introduction of proper initial (equilibrium /o( r ) p) an d initial perturbation /(r, p, 0)) and 
boundary conditions for the quantum distribution function (Wigner function) /(r, p, t) 45 1. 
Generally, these conditions differ from those for the classical distribution function, since the 
Wigner function is not a positively defined quantity and hence is not equivalent to particle 
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density in coordinate- momentum phase space 32|. The equilibrium quantum distribution 
function should take into account the proper quantum statistics (occupation probabilities 
for particle states with given energies) for the considered type of quantum particles (e.g., the 
Fermi-Dirac statistics for electrons, which is not consistent with the approximation of spinless 
particles used in this paper, but still yields correct results, e.g., when modeling electrostatic 
phenomena in electron plasma). We note that the equilibrium distribution function could in 
principle be defined by inclusion into the kinetic model of the collision integral for quantum 
particles St[f] that accounts for both Coulomb and exchange interaction of particles, thus 
automatically ensuring Fermi or Bose statistics for fermions or bosons, respectively. With 
the proper St[f] included in the model, the equilibrium distribution function /o( r >p) °f 
quantum plasma would be defined self-consistently from the equation 



P df 

^- + e 

m or 



F , PxB 



mc 



^-I q0 (v,p) = St[f }. (20) 



However, inclusion of the proper quantum collision integral St[f] greatly complicates the 
kinetic model. Luckily, for many quantum plasmas (e.g., conduction electrons in metals) the 
characteristic collision frequency of particles is many orders of magnitude smaller than the 



28 



so that collisions are 



characteristic frequency of plasma dynamics (plasma frequency) 
not important for plasma dynamics (but are important for its equilibrium). In such cases, 
one can still use the collisionless kinetic model for the plasma dynamical response (thus 
avoiding the complications of introducing the collision integral), while postulating (rather 
than obtaining self-consistently from Eq. (|20|) ) the equilibrium distribution function f that 
mimics the effect of the missing St[fo] on defining the equilibrium. The postulated fo 
should be constructed from the wave functions of stationary states of plasma particles in the 
"equilibrium" electromagnetic field defined by <f)Q and An, accounting for the occupational 
probabilities for the corresponding states; see, e.g., Ref. 



2. "Quantum interference integral" 

Introduction of I q (r, p, t) that acts like a source term in the quantum-modified "Vlasov" 
kinetic equation for the quantum distribution function, redistributing /(r, p, t) in coordinate- 
momentum phase space, while preserving the spatial number density of particles according 
to (I18p . Note that the dimensionality of the integral in I q (r, p, t) is significantly reduced 
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by using symmetries of a chosen ad hoc problem, from the general six- dimensional integral 
down to a two-dimensional integral for one-dimensional problems, and to a four- dimensional 
integral for two-dimensional problems. Thus the evaluation of I q {r, p, t) is not expected 
to be prohibitively computationally expensive, especially for one-dimensional problems. In 
regimes when quantum effects are small but finite (i.e., when h is small compared with all 
plasma parameter combinations with dimensionality of %\ e.g., em^n -1 / 6 ), approximation 
(fT9]) for I q (r, p, t) should be used instead, which can be further simplified by using symmetries 
of the problem (e.g., in one- dimensional case the integration in ([19]) can be done analytically 
to the end). 

These modifications are convenient since they do not affect the main parts of the classical 
Vlasov code: advection and acceleration of f(r, p) in coordinate-momentum phase space, 
evaluation of charge and current densities, and solution of equations ©-([7]) for the scalar and 
vector potentials of the self-consistent electromagnetic field, - all remain unchanged, owing to 
the formal similarity of the quantum kinetic equation ( fl6l) to the classical Vlasov-Boltzmann 
equation. Thus the suggested modification of an existing classical electromagnetic Vlasov 
code 42h44|. instead of developing an electromagnetic analogue of Wigner-Poisson code 46] 
that would numerically solve Eq. @ with computationally prohibitive integral on the right 
hand side, offers, in our view, a much quicker and more straightforward path to simulating 
linear and nonlinear electrostatic and electromagnetic phenomena in nonrelativistic quantum 
plasmas. 



B. Linear response of uniform isotropic quantum plasma 

To obtain the linear dielectric response of a quantum plasma, we consider a small per- 
turbation f(r,p,t) to equilibrium characterized by quantum distribution function fo(r,p) 
formed in the presence of stationary electromagnetic field characterized by "equilibrium" 
potentials 0o an d A . (We note again that fo can not be obtained self-consistently from the 
kinetic model (jl~6l) . but rather needs to be postulated (constructed from stationary states of 
plasma particles in the "equilibrium" field defined by 4> and A ), as discussed in Sec. lIII All ) 

In case of a uniform isotropic equilibrium, we have (fio = A = 0, /o(r,p) = fo(p). 
Substituting /(r, p, t) = fo(p) + /(r,p,t) with |/| f into ( 1T6]) and linearizing, one 
obtains, assuming V/o = (uniform equilibrium) and V ■ A = (Coulomb gauge), the 
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following linear equation for /(r, p, t): 



df , P df 

"aT H a - + e 

ot m ot 



„ , 1 OA p x V x A 

-V0-- -~r + - 

c OT mc 



9/o ie 



dp (2ir] 



dAdfe iA -«- p) /o(0 



x ■ ^ - 4, «)]) - 1 V ' T^ 1 " ^ ' A <" «» I > ■ <«> 



Fourier transforming (}2~Tj) . integrating over A and £, and solving for /, we obtain fit, r, p) = 
(2-7r)~ 4 JJ dujdkfu^p) exp[— i(ut — k ■ r)], with 

A,,k(p) = (-e^, k + — (P ■ A^)) Dp ' k[/ °) + — (P " A., k ) ^° , (22) 
V mc / a; — k • p/m mc ae 

where Wj k and A Ui k are the Fourier transforms of 0(r, £) and A(r, i), respectively, e = p 2 /2m 

is kinetic energy of plasma particles, and the "momentum difference" operator -D p ,k is defined 

as £>p, k [/o(p)] = L/o(p + ftk/2) -/ (p-7ik/2)]/7i. Finally, substituting the obtained A,, k (p) 

into the Fourier transformed Eq. (J7J), after some straightforward manipulations one obtains 

the linear dielectric permittivity tensor £ij(cu, k) of uniform isotropic quantum plasma in the 

form 

£i >, k) = s\u, k)^ + e*(u, k) - ^ , 
with the longitudinal and transverse permittivities given by 

e(w ' k) = i+ ^y ^- k .p/ m ' (23) 

k) = 1 - 4 + ^ / ^P pi"^-, (24) 
m z oo z J 00 — k ■ p/m 



where is the absolute value of the component of p perpendicular to k. 

Note that in (123l) - (124l) the integration over p\\ (a component of p along k) is undefined 
due to a singularity (simple pole) at p\\ = mu/k, for real uj,k. This difficulty is avoided 
by requiring that all perturbations obey the causality principle, i.e., appear at some initial 
time to- F° r such perturbations, their temporal Fourier transforms become one-sided (e.g., 
(fiui = J °° <t>{t') exp(iu;t')dt' , with t' — t — to), and are defined for complex uj with Im(u;) > 0. 
To extend their definition to real u, one needs to perform their analytical continuation 
onto the real axis lm(u) = of complex u> plane, which can be done by taking the limit 
Im(oj) — > 0+ in the corresponding functions of complex u. This procedure, known as 
Landau's rule, defines how the contour of integration over p\\ in (|23|) - (|24|) should be deformed 
to avoid the singularity at p» = mu/k, and is equivalent to replacing to with Re(a;) + io 
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and taking the limit o — > 0+, followed by integrating along real pn. The resulting dielectric 
permittivity can be complex, its imaginary part leading to Landau damping 
The obtained linear longitudinal and transverse dielectric permittivities 



33]. 



match 

with those previously obtained in the literature for uniform isotropic weakly coupled quan- 



tum plasma of spinless particles 



20|, and define, e.g., dispersion properties of small- 



amplitude electrostatic and electromagnetic oscillations in such plasmas. 



C. Applicability 

The quantum kinetic equation (1TB"]) is obtained for weakly coupled nonrelativistic quantum 
plasma. The assumption of weak coupling is justified when the coupling parameter T = 
f^int/ekin is small, which in case of plasma consisting of degenerate electrons with density n 
and Fermi energy ep = ( y h 2 /2m)(3n 2 n) 2 / 3 implies T ~ e 2 n 1 ' 3 /ep <C 1, which is equivalent to 

e 2 /hv F < 1, (25) 



where vp = \j2epjm is the electron Fermi velocity. (Note that, even if the spin of plasma 
electrons is taken into account, the corresponding effect of electron correlation due to their 



exchange interaction is negligibly small when ( 125]) is satisfied |20|.) On the other hand, the 
assumption of nonrelativistic plasma, v <C c, where v is the characteristic velocity of plasma 
particles, and c is the speed of light, implies in the case of degenerate electron plasma (for 
which v ~ Vf) that 

e 2 /hv F > e 2 /hc « 1/137. (26) 

If the condition (1231) is violated (e.g., for electrons in metals), the effect of particle corre- 
lations (collisions and exchange interactions) becomes important, and such plasmas should 
be treated as a quantum liquid; however, the present model (fTBT) may still yield qualita- 



tively correct results for such plasmas [20|. In plasmas where condition (126]) is violated 
(e.g., in dense plasmas subject to ultrarelativistic laser radiation, or in the matter of white 
dwarfs), relativistic effects become important. Both conditions (1231) and ( 123]) may be met 
simultaneously in some semiconductors. 

The model presented here ignores the spin of plasma particles, which makes it formally 
valid only for quantum plasmas consisting of spinless particles, with conditions (1231) and 
(T23]) simultaneously met. On the other hand, our kinetic model is more general than any 

14 



quantum fluid model, as it accounts for essentially kinetic effects which the fluid models 
ignore completely, and thus it should be used in cases where such effects are expected to be 
important. Generalizing our kinetic model to include either (or, ideally, all) of the ignored 



effects (spin ll|, |20|, plasma particle correlations 20], and relativistic effects) would be a 



natural continuation of this work. 



D. On gauge invariance of the model 



The Wigner function /(r, p, t) used in the kinetic model presented here explicitly depends 



30 



40|. Indeed, applying the 



on the potentials and A, and thus it is not gauge-invariant 
gauge transformation of potentials A' = A + Vg(r,£), <$' = <p — (l/c)dg(r,t)/dt along with 
the simultaneous unitary transformation ip' = ip exp[(ie/hc)g(r, t)} of particle wave functions 
(introduced in order to preserve the form of Schrodinger equation with gauge-transformed 



potentials 
/'(r,P,t) 



4l|), we have for the transformed Wigner function: 



(2vr) : 



ciAexp 



le 



-(\.Vg-D rX [g\ 



-iA-[p+(e/c)A] ( _ 



p(r -h\/2,r + h\/2). 

(27) 

— * ^ 

In the classical limit, A ■ Vg — D r ^[g] vanishes, leading to /' = /, i.e., the Wigner function 
becomes the gauge-invariant classical distribution function, as expected. However, in general 
/'(r, p, t) 7^ /(r, p, t), i.e., the quantum distribution function /(r, p, t) is not gauge- invariant. 

The fact that /(r, p, t) is gauge-dependent presents a problem of defining high-order (3rd 
and higher) moments of / in a gauge-invariant way 30j . This problem may be important in 
the hydrodynamic model of quantum plasma, where the evolution of lower-order moments 
of / relies on (and is affected by) the evolution of higher-order moments, which thus all 
have to be defined from /(r, p, t) in a gauge-invariant way as the chain of hydrodynamic 
equations is derived from the kinetic model {(], [joj]. However, in the kinetic description of 
plasma electrodynamics developed here, the problem of gauge-dependent /(r, p, t) does not 
present itself. Indeed, in order to have a complete description of self-consistent electromag- 
netic field in plasma, only the knowledge of charge and current density of plasma particles 
is required, as only these two quantities appear in the Maxwell's equations for the self- 
consistent electromagnetic field. In the kinetic description employing the gauge- dependent 
Wigner function /(r, p, t), electron charge and current densities, defined by Eqs (E])-® 
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as the zeroth- and first-order moments of /, are gauge-invariant, even though / itself is 



not 



30| ; hence the plasma response to the electromagnetic field, defined from the gauge- 



dependent Wigner function, is gauge-invariant. Therefore, the kinetic equation (|T6"|) coupled 
with Maxwell's equations for <f) and A (written for a chosen gauge), although using the 
gauge-dependent Wigner function /(r, p, t), describes electrodynamics of quantum plasma 
in a gauge-invariant way. (We stress that Eqs (fl6]) - (IT?]) are derived without assuming any 
gauge fixing, and thus can be used in their present form for any gauge.) 

In those cases where knowledge of third- or higher-order moments of /(r, p, t) is required 
from a kinetic simulation, an alternative kinetic model employing a specially constructed 
gauge-invariant Wigner function may be of advantage 30j. Otherwise, for most practical 
purposes, the kinetic model developed here provides an adequate description of quantum 
plasma electrodynamics, subject to the applicability limits discussed in Sec. IIII CI 



IV. CONCLUSIONS 



In this paper, a kinetic theory of nonrelativistic quantum plasma with electromagnetic 
interaction of its particles is developed in the framework of the Hartree's mean-field approx- 
imation. The obtained quantum kinetic equation for the one-particle quantum distribution 
function (Wigner function) is cast in a form of the classical Vlasov equation with the addi- 
tional quantum source term accounting for the quantum interference of overlapping plasma 
particle wave functions. This form of quantum kinetic equation suggests a straightforward 
modification of existing classical electromagnetic Vlasov-Maxwell codes, without affecting 
the core parts of their solvers, which opens a relatively quick path to simulating electro- 
magnetic (as well as electrostatic) phenomena in various quantum plasma structures. Such 
modified Vlasov codes could also provide a useful basis for theoretical studies of quantum 
plasmas, as quantum and classical effects can be easily separated in such codes. 
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